Pokhilko et al. BMC Systems Biology 201 3, 7:23 
http://www.biomedcentral.eom/1752-0509/7/23 



Systems Biology 



RESEARCH ARTICLE Open Access 



Modelling the widespread effects of TOC1 
signalling on the plant circadian clock and its 
outputs 

Alexandra Pokhilko 1,4 , Paloma Mas 2 and Andrew J Millar 1,3 " 



Abstract 

Background: 24-hour biological clocks are intimately connected to the cellular signalling network, which 
complicates the analysis of clock mechanisms. The transcriptional regulator TOC1 (TIMING OF CAB EXPRESSION 1) is 
a founding component of the gene circuit in the plant circadian clock. Recent results show that TOC1 suppresses 
transcription of multiple target genes within the clock circuit, far beyond its previously-described regulation of the 
morning transcription factors LHY (LATE ELONGATED HYPOCOTYL) and CCA1 (CIRCADIAN CLOCK ASSOCIATED 1). It 
is unclear how this pervasive effect of TOC1 affects the dynamics of the clock and its outputs. TOC1 also appears to 
function in a nested feedback loop that includes signalling by the plant hormone Abscisic Acid (ABA), which is 
upregulated by abiotic stresses, such as drought. ABA treatments both alter TOC1 levels and affect the clock's 
timing behaviour. Conversely, the clock rhythmically modulates physiological processes induced by ABA, such as 
the closing of stomata in the leaf epidermis. In order to understand the dynamics of the clock and its outputs 
under changing environmental conditions, the reciprocal interactions between the clock and other signalling 
pathways must be integrated. 

Results: We extended the mathematical model of the plant clock gene circuit by incorporating the repression of 
multiple clock genes by TOC1, observed experimentally. The revised model more accurately matches the data on 
the clock's molecular profiles and timing behaviour, explaining the clock's responses in TOO over-expression and 
tod mutant plants. A simplified representation of ABA signalling allowed us to investigate the interactions of ABA 
and circadian pathways. Increased ABA levels lengthen the free-running period of the clock, consistent with the 
experimental data. Adding stomatal closure to the model, as a key ABA- and clock-regulated downstream process 
allowed to describe TOC1 effects on the rhythmic gating of stomatal closure. 

Conclusions: The integrated model of the circadian clock circuit and ABA-regulated environmental sensing allowed 
us to explain multiple experimental observations on the timing and stomatal responses to genetic and 
environmental perturbations. These results crystallise a new role of TOC1 as an environmental sensor, which both 
affects the pace of the central oscillator and modulates the kinetics of downstream processes. 
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Background 

Circadian clocks allow most eukaryotes and some 
prokaryotes to anticipate the environmental day/night 
cycle, through rhythmic modulation of multiple physio- 
logical processes [1]. Circadian clocks are characterised 
by free running rhythms with a period of -24 h even in 
the absence of any environmental cues, such as in 
constant light conditions. Their timing is synchronised 
with the environmental day-night cycle principally by 
responses to light and temperature, but clock circuits 
also respond to many additional stimuli. The signalling 
pathways from such stimuli are often found to be rhyth- 
mically controlled by the clock, forming nested feedback 
loops that modify the circadian oscillator. The gene 
circuit of the higher plant clock consists of interlocked 
morning and evening loops operating through negative 
feedback mechanisms [2]. TOC1 is an important 
component of the plant clock since its mutation or 
overexpression dramatically change the properties of 
the central oscillator [3]. 

Although the importance of TOC1 in the plant circa- 
dian clock has long being recognized, the precise role of 
TOC1 in the clock has been a matter of debate [2-7]. 
Based on indirect genetic evidence it was first proposed 
that TOC1 activated the expression of the key transcrip- 
tion factors LHY and CCA1 [4], although several obser- 
vations on clock mutants remained paradoxical. Recent 
experimental and computational data [2,3,7] have dem- 
onstrated that TOC1 functions as a repressor rather 
than an activator. The studies also show that TOC1 
repression is not limited to CCA1 and LHY but to 
nearly all of the genes at the core of the oscillator [3]. 
The observed reduction in the expression of multiple 
clock genes by up-regulation of TOC1 could not be 
described by our previous model, where only LHY and 
CCA1 were affected by TOC1 [2,3]. Here, we extend 
our previous model by including the repression of mul- 
tiple clock genes by TOC1 and explore TOC1 effect on 
the clock. In particular, we analysed clock dynamics in 
various TOCl-misexpressed lines and addressed the 
long-standing paradox of simultaneous reduction of 
LHY and CCA1 amplitudes in both tocl mutant and 
TOCl-overexpressed (TOCl-ox) plants [7,8]. 

The recently discovered pervasive negative effect of 
TOC1 on gene expression [3,7], the overlap between 
microarrays of the TOC1- and ABA-regulated genes and 
the gated induction of TOC1 by the stress hormone 
abscisic acid ABA [9] suggests that TOC1 might integrate 
environmental and circadian information to regulate 
downstream physiology at specific time of the day [9]. 
Here we combined the clock model with key processes 
regulated by ABA and explored mutual interactions 
between clock and ABA pathways. The model describes 
the regulation of stomatal aperture, one of the most 



important outputs of ABA signalling, which is gated by 
the clock [9]. Being a sensor of water availability in plants, 
ABA level increases under dry conditions, which acceler- 
ates closing of stomata in the afternoon [9-11]. In agree- 
ment with the microarray studies, TOC1 mis-expressing 
plants were shown to have altered ABA-dependent sto- 
mata closure and affected responses to drought conditions 
[9]. To model the observed gating of stomatal opening by 
the clock, we included the interaction of ABA with the 
ABAR protein (also known as the H subunit of Mg- 
chelatase, CHLH or GUN5) - a key component of the 
ABA signalling pathway, arguably postulated as ABA re- 
ceptor [12]. Besides being one of the downstream targets 
of TOC1, ABAR was shown to be a key component in the 
gating of ABA signalling by the clock [9,13]. We used the 
model to explore the dynamics of the opening of stomatal 
pore under various genetic and environmental perturba- 
tions, such as mutations of clock genes and by changing 
the light and water conditions, and showed that the model 
provides good match to existing experimental data. Our 
results show how the dual roles of TOC1 are integrated, 
combining timing and environmental information to 
modulate downstream physiology. 

Results and discussion 

Model description 

The latest plant clock model [2], herein called P2011, 
was revised by adding the recently discovered negative 
regulation of the oscillator genes by TOC1 (blue lines 
on Figure 1A) as described in the section 1 below. 
Additionally, we have included a simplified version of 
the main steps leading to the induction of TOC1 by 
ABA and the regulation of stomata aperture (Figure IB), 
described in the section 2 below. 

1. Extending the clock model by inhibition of gene 
expression by TOC1 

The latest model of the Arabidopsis oscillator [2] 
was extended based on our recent data on the 
inhibition of expression of multiple clock genes 
by TOC1 and induction of TO CI through ABA 
signalling. The model is described by a system of 
32 ODEs. The basic structure of the plant clock 
was kept from our previous P2011 model [2]. 
Briefly, the key components of the clock are the 
LHY and CCA1 morning genes and the evening 
complex (EC) genes. I/fFand CCA1 are 
similarly expressed around dawn and act 
together to regulate multiple targets. They are 
described by single variable LHY/CCA1 as 
before. EC (EVENING COMPLEX) genes LUX, 
ELF3 (EARLY FLOWERING 3) and ELF4 are 
expressed around dusk and form the EC protein 
complex, which suppresses expression of multiple 
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The previous P201 1 model was extended by adding: 

1) TOC1 inhibition terms to the equations for PRR9, PRR7, Gl, LUX and 
ELF4 transcription. 

2) Induction of TOC1 transcription by ABA signalling through SnRK2 

3) The key components of ABA signalling that are related to the clock: ABAR, 
PP2C, SnRK2 and stomatal aperture 

Figure 1 Extended model of the plant circadian clock, including the multiple targets of TOC1 and the interactions with the ABA 
signalling pathway. A: The clock model was extended from the previous one (P201 1) by including the negative transcriptional regulations of 
the core clock genes by T0C1 (blue lines) and the up-regulation of TOO expression by ABA signalling, as shown in panel B. Elements of the 
morning and evening loops are shown in yellow and grey, respectively. Proteins are shown only for EC, ZTL and C0P1 for simplicity. Post- 
translational regulation of T0C1 and the EC by Gl, ZTL and C0P1 is indicated by the dotted red lines. B: Simplified model of ABA signalling, gated 
by the clock through regulation of ABAR expression. See section "Model description" for more details. Transcriptional and protein regulation are 
shown by solid and dotted lines, respectively. C: Short description of the newly-introduced connections and elements of the model compared 
to P2011. 



genes at night (Figure 1; [2]). LHY/CCA1 
expression is suppressed by the PRR proteins in 
a morning loop, while the EC is negatively 
regulated by the ubiquitin E3 ligase COP1 
(CONSTITUTIVE PHOTOMORPHOGENIC 1), 
which targets ELF3 protein to degradation by 
proteasome and by Gl protein, in the evening 
loop. The morning and evening loops are further 
interlocked through the suppression of evening 
gene expression by LHY/CCA1 and suppression 
of PRR9 and TOC1 (PRR1) expression by the EC. 
The model also includes the F box protein ZTL 
(ZEITLUPE), which negatively regulates the level 
of TOC1 protein. Light entrains the clock 
similarly to the P2011 model through several 
mechanisms, which are supported by experimental 
data [2]. The main mechanisms include acute 
activation of LHY/CCA1, PRR9 and GL transcription 
(eqs. 1, 5, 27) immediately after dawn; stabilization 
of PRR proteins in presence of light (eqs. 6, 8, 10) 
and light-dependent regulation of the EC by COP1 
and Gl proteins (eqs. 24, 21, 22). 



The TOCl-related reactions were substantially 
extended. Firstly, we added multiple reactions of 
inhibition of clock gene expression by TOC1. 
This includes direct inhibition of expression of 
LHY/CCA1 as in the P2011 model (eq. 1), and 
also PRR9, PRR7, PRRS, LUX, ELF4 and GL (eqs. 
5, 7, 9, 18, 14, 27), which is based on our data 
[3] (Figure 1A). The details of TOC1 interactions 
with other regulators (e.g. LHY/CCA1, EC, 
protein P) to modulate gene expression are 
largely unknown, so we assumed that TOC1 acts 
as a non-competitive inhibitor of gene expression 
Secondly, we added physiologically relevant 
activation of TOC1 expression by the ABA 
signalling pathway (Figure IB, [9]), described in 
section 2 below. This provided an additional 
level of input to the clock through ABA, which 
is induced by stress. Thirdly, we added 
regulation of stomata aperture as a clock output, 
which is directly related to both ABA and TOC1 
signalling (Figure IB), as presented in section 2 
below. 
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The model is described by the following system of 
ODEs: 
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Where c™ and stand for dimensionless 
concentrations of mRNA and protein, respectively. 
The time unit is an hour. Index "i" labels the 
molecular components as follows: 

ABAR ABAR mRNA 

AR ABA complex with ABAR protein 



COPlc COP1 cytoplasmic protein 

COP Id COP1 nuclear protein, day 

COP In COP1 nuclear protein, night 

E3 ELF3 mRNA 

E3c ELF3 cytoplasmic protein 

E3n ELF3 nuclear protein 

E4 ELF4 mRNA and protein 

E34 nuclear protein complex ELF3-ELF4 

EC nuclear protein complex ELF4-ELF3-LUX 

EGc cytoplasmic protein complex ELF3-GI 

EGn nuclear protein complex ELF3-GI 

G GI mRNA 

Gc GI cytoplasmic protein 
Gn GI nuclear protein 

Gn_tot total amount of GI protein in nucleus 

L LHY/CCA1 mRNA and protein 

Lmod LHY/CCA1 modified protein 

Ltot total amount of LHY/CCA1 protein 

LUX LUX mRNA and protein 

PS PRR5 mRNA and protein 

P7 PRR7 mRNA and protein 

P9 PRR9 mRNA and protein 

PP2C PP2C active protein 

SnRK2 SnRK2 active protein 

T TOC1 mRNA and protein 

ZG cytoplasmic protein complex GI-ZTL 

ZTL ZTL protein 

Symbol "s" corresponds to stomata aperture in 
relative units (the maximum of s is equal to 1). A 0 is 
the total ABA level in relative units. ABA level for 
most of calculations was set to 1, unless it is stated 
otherwise (ABA 0 =1). 

The quasi-steady state approximation for the 
nuclear complexes ELF3-GI, ELF3-ELF4 and 
nuclear GI protein, which reduces the number of 
equations, was taken from the P2011 model. 
COP1 nuclear activity was described by two 
forms of COP1 ([2], eqs. 20, 21). The parameters 
rij represent the rate constants of transcription 
and stomata opening, m ; are the rate constants 
of degradation; pj are constants of translation, 
protein modification and protein complex 
formation; gj are Michaelis-Menten constants and 
a, b, c, d, e,f, g, h, i, j are Hill coefficients; are 
the rate constants of acute (P-dependent) light 
activation of transcription. The acute light 
response in activation of PRR9, LHY/CCA1, GI, 
transcription (eqs. 5, 1, 27) was modelled using a 
light-sensitive activator - protein P (c P ), which is 
accumulated in darkness and was degraded in 
light similarly to the P2011 model (eq. 4). L= 1 
when light is present, 0 otherwise; D = 1-L. TheL 
(t) function was used to simulate smooth 
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transitions between L and D analogous to the 
P2011 model: 

L(t) = 0.5-((l + tanh(0 - 24^?oor(f/24) - dawn)/T)) 
-(1 + tanh(0 - 24-floor(t/24) - dusk)/T)) 
+ (1 + tanh((£ - 24;/ZoorO/24) - 24)/ r))) 

Where dawn and are the phases of dawn and 
dusk (normally dawn = 0); T is the duration of 
twilight (we used T = 0.05 h); tanh and floor - 
standard functions of hyperbolic tangent and 
rounding operation. This representation of periodic 
input corresponds to the Input Signal Step 
Function, which now has convenient software 
support in SBML [14]. Mutations of the clock genes 
were simulated by decreasing the rate of 
transcription of the corresponding gene to zero. 
TOCl-ox plants were modelled by adding the 
background activation of TOC1 translation 
(parameter nl6 in eq. 12; nl6 = 0 for wild type). 
The equations for COP1 E3 ligase activities 
(eqs. 20, 21), which were developed in P2011, are 
independent of the clock and their parameters were 
constrained from the data on the kinetics of COP1 
substrates HFR1 and HY5 proteins upon dark/light 
transitions similarly to P2011 [15]. The equations for 
HY5 and HFR1 proteins, which were used only for 
the optimization of COP1 parameters, are taken 
from P2011: 

dc H Ys/dt =p 2 2 - m^C H YsCC0Pld 
- m 2 5C H Y5C C OPln 

dCHFRl/dt =P30 - rn 2 8CHFRlCC0Pln 

2. Modelling of ABA signalling and its relation with 
the clock 

To model the mutual connections between the clock 
and ABA signalling, we extended the clock model by 
including ABA-induced reactions. 
The first group of reactions of ABA signalling is 
related with the inactivation of protein phosphatase 
PP2C (protein phosphatase of 2C type) by the ABA 
complex ([11]; Figure IB). Based on the data 
showing relatively small diurnal changes in ABA 
concentration in plants [10,16], we assumed that 
ABA level is constant in our model and treated it as 
a parameter ABA 0 . Next we included the potential 
ABA receptor ABAR, which is a necessary 
component of the regulation of both TOC1 
expression and stomatal aperture by ABA [9]. ABAR 
expression has a strong circadian pattern, peaking in 
the morning [9]. Our analysis of the ABAR 
promoter in Arabidopsis revealed 2 CCA1 binding 
sites (CBS) AAATCT [17], in addition to the known 



TOC1 -binding sites [3,9], suggesting a possible 
mechanism for the regulation of ABAR transcription 
by the clock. This was described in our model 
through activation of ABAR expression by LHY/ 
CCA1 and its inhibition by TOC1 ([3,17], eq. 29). 
To simplify the model and reduce the number of 
equations, we used a quasi-steady state 
approximation for the ABA complex with ABAR, 
assuming that ABAR protein level follows the 
cognate mRNA. PP2C activity was described 
through its inhibition by this ABA complex (eq. 30). 
The second important group of reactions of ABA 
signalling is related with the activation of SNF1 
(Sucrose-Nonfermenting Kinasel)-related protein 
kinase SnRK2 by ABA [11]. This activation happens 
through double-negative interactions: the ABA 
complex inactivates PP2C and PP2C inactivates 
SnRK2. So we included in the model an equation for 
SnRK2 activity, which is inhibited by PP2C (eq. 31). 
After activation by ABA, SnRK2 mediates multiple 
reactions of the ABA pathway, including induction 
of transcription of multiple target genes which 
possess ABRE elements (ABA-responsive elements) 
in their promoters [11]. As TOC1 was reported to 
be induced by ABA [9], we analysed the TOC1 
promoter and identified 2 ABRE elements. This 
suggested that a mechanism of TOC1 induction by 
ABA is similar to other genes and mediated by 
SnRK2 (Figure IB). Therefore we included in the 
model the activation of TOC1 expression by SnRK2, 
assuming that this activation competes with 
inhibition of TOC1 by LHY/CCA1 protein ([9], 
eq. 11). An analogous mechanism of ABA signalling 
through the ABA/PP2C/SnRK2 pathway is involved 
in the acceleration of stomatal closing by ABA [11], 
which we included into the model. The equation for 
stomatal aperture describes stomatal opening in a 
phenomenological way, which allowed us 
qualitatively compare the model behaviour with 
physiological data. Additionally to the inhibition of 
stomatal opening by ABA signalling through active 
SnRK2, we included the known activation of stomatal 
opening in the presence of light ([18], eq. 32). 
So finally the clock model was extended by 
including eqs. (29)-(32) for ABAR, PP2C, SnRK2 and 
stomata aperture (variable "s"). Thus we connected 
the clock with a simple version of ABA signalling, 
which enabled us to simulate TOC1 and stomatal 
regulation by ABA. The effective parameters of the 
ABA-related reactions were optimized together with 
other parameters the clock model against multiple 
perturbations as described below. 
Additional file 1: Figure SI shows simulated diurnal 
kinetics of the main components of ABA signalling 
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under the optimal parameters, which are presented 
in Additional file 2: Table SI. The peak of ABAR 
transcription soon after dawn corresponds to the 
publicly available microarray data (http: //diurnal. 
mocklerlab.org/), which agrees with RT-Q-PCR 
analysis [9] . The accumulation of ABA- ABAR 
complex in the day time leads to inactivation of 
PP2C and activation of SnRK2, with the peak of 
SnRK2 activity in the afternoon (Additional file 1: 
Figure SI). SnRK2 then induces TOC1 expression 
and stimulates the closing of stomata as described 
above. 

3. Solving the system of ODEs 

The extension of P2011 by adding newly discovered 
processes of the inhibition of the PRR9, PRR7, PRR5, 
LUX, ELF4 and GI clock genes by TOC1 (Figure 1A) 
resulted in the addition of 7 new parameters to 
P2011. Our model was also extended compared to 
P2011 by adding reactions describing ABA-related 
processes (Figure IB), which introduced 14 more 
new parameters. Finally, the model describes the 
dynamics of stomatal opening, which was absent 
from P2011, resulting in 5 more new parameters. 
Therefore, the higher level of biological complexity 
of the model and its extension to hormonal 
signalling through ABA inevitably increased the 
parameter space, adding 26 parameters compared to 
P2011. This allowed our model to describe 
experimental data that were not described by P2011, 
such as the inhibition of multiple clock genes by 
TOC1, stimulation of TOC1 expression by ABA 
signalling, changing the clock period by ABA and 
the dynamics of stomata presented in the Results 
below. In addition, the new structure of the clock, 
which now incorporates ABA signalling to the clock 
through TOC1 (Figure 1A) allowed us to improve 
the description of the data on the clock kinetics 
compared to P2011, as shown in the Results below. 
Parameter values were either constrained based on 
experimental data or fitted to multiple time-series 
data sets, similarly to the P2011 model. 49 out of the 
total 133 parameters were constrained based on the 
available experimental data. Most of these 
parameters have the same values as in P2011 
because they were constrained by the same data. 
Other parameters were fitted to multiple datasets 
similarly to P2011, but their values may be different 
from P2011 because the structures of the models are 
different. The optimal set of parameters is presented 
in Additional file 2: Table SI. The new model, under 
the optimal parameter values, retains most of its 
properties from P2011, keeping the good fit to data 
from various conditions, such as diurnal cycles of 
the wild type and Ihy/ccal, Ihy/ccal/gi and elf3 



mutants, as well as continuous light or darkness in 
the wild type and multiple mutants. For example, 
the period values for the simulated free-running 
clock in wild type and mutant plants were 24.5 h for 
wild type under constant light conditions, 26.6 h for 
wild type under constant darkness and 17.7 h, 
21.4 h, 30.1 h, 30.7 h, 21.2 h for Ihy/ccal, tocl, ztl, 
prr79 and gi mutants in constant light, which match 
the experimental observations [8,19-22]. 
Additionally Ihy/ccal/gi, elJS, elf4 and lux mutants 
were arrhythmic in constant light, in agreement 
with the data [8,23-25]. We also tested the 
sensitivity of the model to the variations of the 
newly introduced parameters of the inhibition of 
target genes by TOC1 and ABA signalling. 
Additional file 1: Figure S2 shows the relative 
changes in the amplitude of LHY expression and the 
clock period in constant light conditions under 10% 
changes of each parameter. We observed less than 
3% changes in the period and less than 12% change 
in the amplitude of LHY mRNA under variations of 
these parameters, which shows that the extended 
model retains its robustness to parameter 
perturbations from P2011 (Additional file 1: Figure 
S2). Thus, the optimal set of parameters 
demonstrated good correspondence to multiple 
datasets and provides robust behaviour in the model. 
However, given the significant complexity of the 
model, which is related to the complexity of the 
biological processes involved, and the number of 
unknown parameters, we cannot exclude the 
existence of other parameter sets that might equally 
describe the data. The detailed study of this large 
parameter space lies outside the scope of this paper. 
The model with the current parameter set matches 
well enough to the data to make it a useful 
instrument for the exploration of the possible 
mechanisms of the observed biological phenomenon. 
ABA modulates expression of the clock genes 
through the induction of TOC1 expression followed 
by suppression of TOC1 target genes. Additional file 
1: Figure S3 shows the dependence of the peak 
mRNA values of all clock genes on the ABA levels, 
demonstrating that ELF4, LUX, GI, PRR7 and PRRS 
are the most sensitive to ABA in our model. 

Modulation of circadian dynamics via clock gene repression 
by TOC1 

The extension of our previous P2011 model with the par- 
allel suppressive effects of TOC1 on the oscillator genes 
(Figure 1A; [3]) allowed us to improve the description of 
data from several clock mutants compared to P2011, with- 
out affecting the good fit of the P2011 model to the data 
of other mutants. Figure 2 shows the profiles of clock gene 
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Figure 2 Clock gene expression profiles in T0C1 mis-expressed lines. A: The simulated profiles of LHY/CCA1 mRNA under constant light 
conditions are shown for wild type (wt), tod and TOG -ox lines by black solid, dashed and grey lines, respectively. TOCl-ox was simulated with 
nl 6 = 0.1 . Simulations were run after entrainment of the clock to 12 L12D conditions to match the data. B: Experimental data for LHY mRNA in 
wild type and tocl plants under constant light conditions are redrawn from [4]; C: Cartoon illustrating the clock gene profiles in TOC1 mis- 
expressed lines relative to the wild type. Coloured lines represent the expression of LHY and its inhibitors PRR9, PRR7, PRR5, TOG, as shown in the 
legend for the wild type. Blue connections on the wild-type panel (left) show the inhibition of PRR9, PRR7, PRR5 by TOCL The release of this 
inhibition in the toc7 mutant (centre) results in the increase in PRR expression, which causes the reduction of the LHY amplitude and, together 
with the absence of TOC1, shortens the clock period in the toc7 mutant. On the contrary, overexpression of TOC1 (right) reduces PRR expression 
and allows LHY levels to rise in the morning, while higher TOC1 delays LHY at the end of the night, lengthening the period in TOCl-ox. 



expression in the simulated TOC1 -mis-expressing plants 
and explains the mechanisms involved. The observed 
decrease in the peak level of LHY/CCA1 expression in the 
tocl mutant suggested that TOC1 was an activator of 
LHY and CCA1 [4]. Although TOC1 is a repressor of 
LHY/CCA1 in the model, removing TOC1 also reduces 
LHY/CCA1 expression (Figure 2A). The simulated tocl 
mutation increases the levels of the remaining LHY/CCA1 
inhibitors, the PRR (PSEUDO-RESPONSE REGULATOR) 
proteins (PRR9, PRR7, PRR5), by relieving the negative 
regulation of PRR transcription by TOC1 (Figure 2C). In 
the model, this results in instant down-regulation of LHY/ 
CCA1 mRNA after switching from a simulated light/dark 
cycle to constant light (LL) conditions. Figure 2A demon- 
strates that the simulated tocl mutant has LHY/CCA1 at 
0.85-fold the peak level in wild type after two days in LL, 
which is close to the experimentally observed 0.82-fold for 
LHY mRNA (Figure 2B, [4,26]). In contrast, the tocl 
mutant simulated in the P2011 model initially contradicts 
the data, as it first shows an increase of LHY/CCA1 levels 
followed by only a slight decrease (0.97-fold after two days 
of LL), and through a more indirect mechanism [2]. 

Over-expression of TOC1 also resulted in a lower level 
of LHY/CCA1 expression, which is due to the direct sup- 
pression of LHY/CCA1 expression by the increased level 
of the repressor, TOC1 (Figure 2A). This effect of TOC1 
was present already in the P2011 model. However, the 



P2011 model failed to describe the longer period of the 
TOCl-overexpressor line {TOCl-ox) under constant light, 
which was observed in the data [27]. The discrepancy in 
the earlier model was caused by increased levels of the 
PRR inhibitors as well as TOC1, which resulted in a re- 
duction of the LHY/CCA1 level. This led to a shortening 
of the period in the simulated TOCl-ox lines in P2011 
rather than a longer period. The current model describes 
the experimentally observed prolongation of the period 
due to delay in LHY/CCA1 in TOCl-ox compared to the 
wild type because repression by TOC1 reduces the level of 
the other PRRs (Figure 2C). The simulated TOCl-ox in 
Figure 2 A has a 1.5 h longer period than the wild type. 
The period difference between the TOCl-ox and wild type 
increases with the increase in the level of TOC1 expres- 
sion in our model, in agreement with the data [21,27]. For 
example, 3-fold increase of the TOC1 overexpression level 
compared to the level shown in Figure 2A results in 5 h 
longer period of the TOCl-ox compared the wild type 
(while it has 1 h shorter period in P2011). The large pro- 
longation of the TOCl-ox period (5 h) corresponds to the 
experimental observations for the TOCl-ox lines with 
high overexpression levels of TOC1 [27]. The oscillations 
are damping in TOCl-ox in our model simulations 
(Figure 2 A) similarly to the data [27], to the point of the 
loss of rhythmicity of the simulated TOCl-ox at very high 
levels of TOC1 -overexpression, which also agrees with the 
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Figure 3 Modulation of clock dynamics by ABA. Diurnal profiles of TOG mRNA (A) under various levels of ABA (ABA level is expressed in 
relative units; the wild-type level is 1). Simulations were run under a 12 L12D diel cycle. B: Free running oscillation of LHY/CCA1 and TOG mRNA 
under constant light conditions in the absence or presence of ABA. The 1.3 h period lengthening in the presence of increased ABA in the model 
simulation is close to the reported 1.5 h lengthening of the period of the clock marker CAB2: LUC under constant light conditions [29]. 



data [21]. The model suggests that the period increase is 
due to the simultaneous up-regulation of TOC1 and 
down-regulation of the other PRR expression in the 
TOC 1 -ox lines, which allows a longer interval of LHY/ 
CCA1 expression than in the wild type. 

Figure 2C shows schematically the mechanisms of 
regulation of LHY/CCA1 expression by the wave of PRR 
inhibitors in lines with different levels of TOC1: In the 
tocl mutant, increased PRR levels reduce LHY/CCA1 
levels, while the absence of TOC1, the last inhibitor, 
shortens the period. In TOC 1 -ox, LHY/CCA1 expression 
is repressed mainly by TOC1 protein, which declines 
when LHY is raising (LHY and TOC1 are anti-phasic) 
and thus allowing oscillations to occur but with longer 
period than in wild type (Figure 2A). Thus the model 
offers an explanation of the low ZJTY and CCA1 levels in 
the tocl mutant [4] that is consistent with the biochem- 
ical data showing that TOC1 is a repressor [3,7]. The 
mechanism depends upon the inhibition of multiple 
clock genes by TOC1, particularly the PRRs (Figure 2C). 

Additionally to affecting the expression of the morning 
genes LHY/CCA1, PRR9, PRR7, PRRS, TOC1 negatively 
regulates the evening genes LUX (LUX ARRHYTHMO), 
ELF4 (EARLY FLOWERING 4) and GI (GIGANTEA) in 
the experimental data [3] and in our model. This is in 
clear contrast with the P2011 model, where TOC1 af- 
fected only LHY/CCA1 expression directly. The new 
model suggests that increased levels of evening gene 
expression result in the arrhythmia of the Ihy/ccal/tocl 
triple mutant under constant light conditions, as opposed 
to the damped oscillation of the Ihy/ccal double mutant. 
This agrees with experimental observations [28]. 

TOC 1 -mediated effects of ABA on circadian timing 

Figure 3A shows simulated diurnal profiles of TOC1 
expression under increasing ABA levels, which might 
correspond to increasing levels of stress, such as drought 



[9,10]. In agreement with the data [9], ABA accelerates 
TOC1 expression in the afternoon. The up-regulation of 
TOC1 under higher ABA levels resulted in the suppres- 
sion of expression of TOC1 target genes in our simula- 
tions (Additional file 1: Figure S3), which agrees with the 
data [3]. The model also suggests that the up-regulation of 
TOC1 by ABA should lengthen the circadian period 
(1.3 h) under free-running conditions (Figure 3B), similar 
to the TOCl-ox line (Figure 2). This agrees with the 
experimental observations, showing periods about 1.5 h 
longer than WT in the expression of the clock output 
gene CAB2 in plants treated with high exogenous ABA 
concentrations [29]. 

Altogether, the extension of the P2011 model including 
the widespread repression of core components by TOC1 
(Figure 1A) and the inclusion of the regulatory function of 
ABA on TOC1 expression allowed us to improve the 
description of multiple datasets on the clock kinetics in 
wild type and mutant plants, and suggest mechanisms to 
explain previously paradoxical data. 

Gating of ABA signalling by the clock 

We next examined the biological relevance of the inter- 
actions between ABA and the circadian clock. To that 
end, we extended the model to the regulation of an 
important physiological process downstream of the 
clock - stomatal aperture. The equation for stomatal 
aperture describes the regulation of stomatal dynamics 
by light and by ABA signals [10,18]. The simulated diel 
kinetics of stomatal aperture are shown on Figure 4, 
together with data from wild-type Arabidopsis plants in 
the same conditions [18]. Light activation results in a 
more open state of stomata during the day than in the 
night, in agreement with the experimental data [18]. 
Stomatal dynamics also have a pronounced diurnal pat- 
tern related to ABA signalling as discussed below. 
Figure 5A shows simulated stomatal dynamics under 
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Figure 4 The kinetics of stomata opening in 8 L16D diel cycles in the model (A) and in the experimental data (B), redrawn from [18]. 

Zero time corresponds to dawn. Solid and dashed arrows on A show the "dusk" and "dawn" anticipations respectively as described in the text. 



different ABA levels. Increased ABA levels correspond 
to drought conditions, which lead to a more closed state 
of stomata, in agreement with the data [10]. The closing 
of stomata by ABA is a well-known physiological re- 
sponse in plants, which is important for the reduction 
of transpiration under drought [10,11]. Additionally to 
the modulation of the average daily size of the stomatal 
pore, ABA regulates the peak and trough phases of 
stomatal aperture. In the model, ABA signal is gated by 
the clock through ABAR transcription in the morning and 
the following gradual increase of SnRK2 (SNFl/Sucrose- 
Nonfermenting Kinase 1 -related protein kinase) activity 
during the day (Figure IB, [11]). The SnRK2 kinase in turn 
regulates multiple downstream processes including the 
dynamics of stomatal aperture [11], as described in the 
Methods. Figure 5A demonstrates that the resulting diur- 
nal timecourse of stomata aperture is characterized by a 
morning peak with some decline towards dusk -"dusk an- 
ticipation", resulting from higher sensitivity of stomata to 
ABA in the afternoon. Similarly, stomata are less sensitive 
to ABA before dawn, which results in an increased stoma- 
tal aperture before dawn - "dawn anticipation". Both 
dawn and dusk anticipations agree with the data [10] and 
are important for closing stomata in the heat of the after- 
noon to reduce a water loss and opening of stomata in the 



cool of the morning to prepare plants for photosynthesis 
after dawn [30]. The modelling here allowed us to 
describe the data on stomatal dynamics and showed that 
circadian regulation of ABAR transcription was consistent 
with the gating of the ABA signalling by the circadian 
clock and the resulting rhythmic stomatal dynamics. 

The circadian regulation resulted in distinctive, free- 
running rhythms of stomatal opening under constant light 
conditions for simulated wild type and mutant plants 
(Figure 5B), similar to those observed experimentally [18]. 
The period of this rhythmic clock output corresponded to 
the one of the clock itself in our simulations, in agreement 
with data [18]. Figure 5B demonstrates this for the short- 
and long-period mutants tocl and prr7/prr9. Additionally 
to the change in period length, the simulated tocl mutant 
had smaller stomatal aperture compared to wild type 
(Figure 5B). To explore the effect of TOC1 on stomata in 
more detail, we compared the dynamics of stomata in wild 
type, tocl mutant and TOCl-ox lines. Figure 6 A demon- 
strates that the simulated stomata were more open in the 
TOCl-ox and more closed in the tocl compared to wild 
type, which fully agree with the experimental observations 
[9]. This effect was due to the direct inhibition of ABAR 
expression by TOCl. Simulation of the model showed 
that this effect was especially pronounced under high 
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Figure 5 Simulated kinetics of stomata opening and its regulation by ABA and the circadian clock. A: Diurnal profiles of stomata aperture 
of wild type plants under 12 L12D conditions for various ABA levels; B: Free running rhythms of stomata aperture under constant light 
conditions, calculated for wild type, tocl and prr7/prr9 mutants. 
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Figure 6 Effect of TOC1 on stomata size. A: Simulated time courses of stomata aperture in wild type, TOCl-ox and tocl mutant under 12 L12D 
cycle for ABA= 1.3. TOCl-ox was simulated at n16 = 0.3. B: Experimental data on stomata conductance in wild type, TOCl-ox and tocl mutant in 
absence (-) and presence (+) of 5 uM of exogenous ABA. Data are redrawn from [9]. 



ABA levels (Figure 6B, [9]). Additionally, the model sug- 
gests that TOC1 regulation will alter the diurnal gating of 
the stomatal response to ABA, because both simulated 
TOCl-ox and tocl mutant lines had less pronounced 
anticipation of dawn and dusk (Figure 6A). 

Conclusions 

The revised clock model explains clock gene profiles of 
mutants with altered TOC1 expression, resulting from 
the pervasive transcriptional repression by TOCl. The 
model was extended to explore the mutual connections 
between TOCl targets and ABA signalling, resulting in 
downregulation of clock gene expression and lengthen- 
ing of the free running period under high ABA condi- 
tions. Moreover, TOCl s participation in the rhythmic 
gating of ABA signalling can explain well-characterised 
rhythms in physiological processes downstream of ABA, 
such as the regulation of stomata aperture. The clock ef- 
fects on stress signalling are significant at the whole- 
plant level, as plant survival under dry conditions was 
reduced in TOCl-ox plants and increased in tocl 
mutants [9]. Moreover, the close overlap between 
TOCl- and ABA-regulated genes in transcriptome pro- 
files [9] suggests that other ABA-mediated stress 
responses are similarly responsive to TOCl. Integrating 
such physiological response pathways, even if they are 
represented in a simplified way, with the clock model 
will allow to expand our knowledge on the temporal 
regulation of plant physiology in its normal context, the 
unavoidable day-night cycle. 

Methods 

The system of ordinary differential equations was 
solved using MATLAB, integrated with the stiff solver 
odel5s (The Math Works UK, Cambridge). The SBML 
version of the model will be available upon publication 



from the Biomodels database [31] (accession number 
BIOMD0000000445) and the Plant Systems Modelling 
portal (http://www.plasmo.ed.ac.uk). A MATLAB version 
of the model is available from the authors upon request. 

Additional files 



Additional file 1: Figure SI. The profiles of the main components of 
ABA signalling, simulated under 12 L12D diel cycles. The daily changes in 
the levels of ABA-ABAR, active PP2C and SnRK2 are shown by black, blue 
and red lines, respectively. Figure S2. The change in the relative 
amplitude and period of LHY mRNA in constant light conditions 
measured under 10% increase and decrease of the model parameters of 
inhibition of target genes by TOC1 and ABA signalling to the clock. 
Figure S3. Relative changes in the peak expression levels of the clock 
genes at varied ABA values. Simulations were done under 12 L12D 
conditions. Expression levels were normalized to the value in absence of 
ABA. 

Additional file 2: Table SI. Optimal set of the model parameter values. 
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